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Abstract: 

We review the status of the leptonic decay constants fx and f}^'^ of the K and K* , respec- 
tively, and the SU(3) breaking quantities ai{K) and af'^(i^'*), the first Gegenbauer-moments 
of the leading- twist distribution amplitudes of K and K* . We obtain new predictions from 
QCD sum rules which are relevant for the calculation of K and K* form factors, for instance 
rpB^K*^ which determines the decay B K*'j, and for QCD factorisation calculations of 
nonleptonic B decays into strange mesons, for instance B — > Kit. 
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1 Introduction 



Hadronic light-cone distribution amplitudes (DAs) of leading twist play an essential role in 
the QCD description of hard exclusive processes. DAs enter the amphtudes of processes to 
which coUinear factorisation theorems apply and were first discussed in the seminal papers 
by Brodsky, Lepage and others [1]. More recently, coUinear factorisation has been shown 
to apply, to leading order in an expansion in l/m^, also to a large class of nonlcptonic B 
decays [2], which has opened a new and exciting area of applications of meson DAs. These 
decays, and in particular their CP asymmetries, are currently being studied at the B factories 
BaBar and Belle and are expected to yield essential information about the pattern of CP 
violation and potential sources of flavour violation beyond the SM. The aim of this letter is 
to provide a reanalysis of SU(3) breaking effects in leading-twist K and K* DAs, using QCD 
sum rules. Our letter is both a sequence to and an extension of previous work reported in 
Refs. [3, 4, 5, 6, 7]. The results are of immediate relevance for all predictions of S — > [K, K*) 
decay processes calculated in QCD factorisation. 

We define two-particle DAs as matrix elements of quark-antiquark gauge-invariant non- 
local operators at light-like separations with z^ = 0. For definiteness we consider dis- 
tributions of mesons with an s quark and a light antiquark q. To leading-twist accuracy, 
the complete set of distributions comprises three DAs (we use the notation z — z^^/j, for 
arbitrary four- vectors z): 

{0\q{z)z-f,[z,0]s{0)\K{q))^ifK{qz) C due-'^^'^''^K{u) , 

Jo 

{0\q{z)z[z,0]s{0)\K*{q,X)) = (e^;.)/! m^,* c^xie-^^^^^^H, 

Jo 

{0\q{z)a,4z,0]s{0)\K*{q,X)) ^ zie';,'^, - ei'\,)fk{fi) f due-^^^^^^ct>i{u), (1) 

^0 



with the Wilson-line 

[2;, 0] = Pexp 



ig I da z'^A^(az) 
Jo 



inserted between quark fields to render the matrix elements gauge-invariant. In the above 
definitions, ei^'' is the polarization vector of a vector meson with polarisation A; there are 
two leading-twist DAs for vector mesons, and (f)j^, corresponding to longitudinal and 
transverse polarisation, respectively. The integration variable u is the (longitudinal) meson 
momentum fraction carried by the quark, u = 1 — u the momentum fraction carried by the 
antiquark. The normalisation constants /^''"""^ are defined by the local limit of Eqs. (1) and 
chosen in such a way that 

»i 

du(j){u) = l (2) 







for all three distributions (f)K, (I>ki ^k- 

The most relevant parameters characterising SU(3) breaking are the decay constants fx 
and /^'", and the first (Gegenbauer) moments of the corresponding leading-twist DAs, a\{K) 
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and tti '"{K*), defined as 

ai{K) = ^ du{2u-l)c})K{u). (3) 

and correspondingly for af'^{K*). ai describes tlie difference of tlie average longitudinal 
momenta of the quark and antiquark in the two-particle Fock-state component of the meson, 
a quantity that vanishes for particles with equal-mass quarks (particles with definite G- 
parity). The decay constants /k and /|^ can be extracted from experiment; has been 
calculated both from lattice [8] and from QCD sum rules, e.g. Ref. [9]. Although the results 
are in mutual agreement, the QCD sum rule calculations are still subject to improvement as 
we shall discuss in this letter. The situation is much less clear with ai. no lattice calculation 
of this quantity has been attempted yet, so essentially all available information on ai comes 
from QCD sum rule calculations. To date, three different types of sum rules have been used 
to calculate oi: 

• sum rules based on the correlation function of two biquark currents of equal chirality, 
so-called diagonal sum rules [3, 4, 6]; 

• sum rules based on the correlation function of two biquark currents of different chirality, 
so-called nondiagonal sum rules [3, 4, 5]; 

• exact operator identities relating ai{K) and a\{K*) to quark-quark-gluon matrix ele- 
ments of the meson, which in turn are calculated from QCD sum rules [7]. 

The basic argument in favour of nondiagonal sum rules is that they are of first order in SU(3)- 
breaking quantities: the perturbative contribution is ~ 0{ms), the leading nonperturbative 
contribution ~ (qq) — (ss), whereas for diagonal correlation functions the corresponding 
contributions are ~ 0{m1) and ms{ss) — mu{uu), respectively. The original calculation of 
Chernyak and Zhitnitsky using a nondiagonal sum rule yielded ai ~ 0.1 [3, 4], but un- 
fortunately suffers from a sign-mistake in the perturbative contribution. This mistake was 
corrected in Ref. [5], which however entails that the two leading contributions come with 
different sign and cancel to a large extent. As pointed out in Ref. [6, 7], the resulting sum 
rules are sensitive to poorly constrained higher-order perturbative and nonperturbative cor- 
rrections and hence numerically unreliable. Alternative sum rules for ai come from diagonal 
correlation functions, a route that was followed, for ai{K), in Ref. [6]. Yet another possibil- 
ity to pin down the elusive ai is offered by exploiting exact operator identities that relate ai 
to quark-quark-gluon matrix-elements which in turn are calculated by QCD sum rules; the 
corresponding results for ai{K) and a\{K*) can be found in Ref. [7]. In this letter we derive 
and analyse diagonal sum rules for the decay constants fx and all three Gegenbauer mo- 
ments ai{K) and a\''^{K*). We plan to come back to the analysis of operator identities and 
the corresponding sum rules for quark-quark-gluon matrix elements in a future publication. 

Out letter is organised as follows: in Sec. 2 we calculate and analyse QCD sum rules for 
Jk, fx^ and ai{K), a\''^{K*). We summarise and conclude in Sec. 3. The appendix contains 
some remarks about the calculation of diagonal sum rules. 
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2 QCD Sum Rules for /^''^^ and a^' ^K^*^) 

QCD sum rules are an established method for the calculation of hadronic matrix elements, 
see Refs. [10, 11] for the original papers and a recent review. The key feature of the method 
is the use of analyticity to relate the local short-distance operator product expansion (OPE) 
of a correlation function of two currents, 

U^i f d'ye^^y{0\TMy)MO)\0) = ^ ROP'^ (4) 

around y = {as opposed to a light-cone expansion aound y'^ = 0, which is appropriate for 
form factor calculations, cf. [12]) valid for = —q^ -C 0, to its dispersion relation in terms 
of hadronic contributions, 

Ti= Tds (5) 

Jo s-q^ -lO 

where p(s) is the spectral density of the correlation function along its physical cut. The OPE 
yields a series of local operators of increasing dimension whose expectation values (0„) in the 
nonperturbative (physical) vacuum are the so-called condensates. In the sum rules analysed 
in this letter, we take into account the condensates and parameters listed in Tab. 1. As for 
the strange quark mass, we would like to recall that with present data there is a hint of a 
discrepancy between unquenched Uf = 2 and rij = 2 + 1 results, the latter ones favouring 
smaller values ms(2 GeV) = (78 ± 10) McV [13]. Awaiting the clarification of this situation, 
we choose to stay with the result from nf = 2 flavours given in Tab. 1. 

The representation of the correlation function in terms of hadronic matrix elements can 
be written as 

p{s) = f6{s-ml) + r^\s), 

where uim is the mass of the lowest-lying state coupling to the currents Ji_2 and p™'^* pa- 
rametrises all contributions to the correlation function apart from the ground state. /, the 
residue of the ground state pole, is the quantity to be determined. A QCD sum rule that 
allows one to do so is obtained by equating the representations (4) and (5) and implementing 
the following (model) assumptions: 



(qq) = (-0.24 ± 0.01)^ GeV'^ 
{qagGq) = (0.8 ± 0.1) GeV^ (qq) 

^^6-2^ = 0.012 GeV^ 


(ss) = (0.8±0.1)(^g) 
{sagGs) = (0.8 ± 0.1){qagGq) 


m,(2 GeV) = (100 ± 20) MeV < — > m,(l GeV) = (137 ± 27) MeV 
Q;,(lGeV) - 0.534 < — > Aqot ° " 384 MeV 



Table 1 : Input parameters for sum rules at the renormalisation scale jj, — 1 GeV. The value 
of nis is obtained from unquenched lattice calculations with Uf — 2 flavours as summarised 
in [13]. 
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• is approximated by the spectral density obtained from the OPE above a certain 
threshold, i.e. p'^™^* p'^^^{s)6{s — Sq) with Sq ~ {itlm + A)^ being the continuum 
threshold, where A ~ O(Aqcd) is an excitation energy to be determined within the 
method. This assumption relies on the validity of semiglobal quark-hadron duality; 

• instead of the weight-functions and l/(s — q^), one uses different weight- 
functions which are optimised to (exponentially) suppress effects of p{s) for large values 
of s and at the same time also suppress higher- dimensional condensates by factori- 
als. This is achieved by Borel transforming the correlation function: Bl/{s — q"^) = 

exp(— s/M^). A window of viable values of the Borel parameter and the 
continuum threshold sq has to be determined within the method itself by looking for 
a maximum region of minimum sensitivity (a plateau) in both and Sq] 

• the OPE of n can be truncated after a few terms. As is well known, this condition is 
fulfilled only for low moments, whereas for higher moments of the DAs in [u — u) the 
nonperturbative terms become dominant. 

After subtraction of the integral over above sq from both sides, the final sum rule reads 

no-"^ ^ ^ £ ds e-^'^' p^^^is) = ^ e-<'^'^\ (6) 

which gives the hadronic quantity / as a function of the Borel parameter and the 
continuum threshold sq (and the condensates and short-distance parameters from the OPE). 
We determine /^''"""^ and ai^''^\K^*^) from the diagonal correlation function 

i j dV"(o|rg(t/)rs(t/)5(o)r[o,^]g(^)|o), (7) 

where z^^ is light-like and the Dirac structures F are given by 

K: r = i75, Kl: F = Kl: T ^ a^^z". 

The calculation with nonlocal operators is very convenient, as it allows one to calculate all 
moments in one go. Specifying for instance to K^,, the sum rule reads 



(8) 



where also n°^^ is expressed as integral over u, which naturally emerges as a Feynman 
parameter in the calculation, and comes with the same weight function exp(w(g2;)). Sum 
rules for fx are obtained as the lowest order in an expansion in qz, those for oi by effectively 
replacing 



-{u-u). 



For f)i^' we find the following sum rules: 

f^^-ml/M^ = SR+, iflfe-^K^^' = SR_, 
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so 



with SR, = -L / + + - fl - 



M2 V 3M2 9 TT y 12M2 ' TT 'V 3 

4 q;^ ms{qq) , 167ra<, w- \ , 167ra^ ,2 , /- \2\ /n^ 
+ 3 V ^ W ^^^^ + W + ^^^^ ) ■ 

This sum rule was aheady given in Ref. [14], apart from the radiative and mass-corrections 
to the quark-condensate contribution which are new. The values of fx and f\ obtained 
from these sum rules are shown in Fig. 1, evaluating all scale-dependent quantities at the 
scale fj, — 1 GeV. The sum rule results agree very well with the experimental values [15] 

fK = (0.160 ± 0.002) GeV, /| = (0.217 ± 0.005) GeV. 

For the situation is slightly more subtle as the correlation function (7) contains contri- 
butions not only from the vector meson K*, but also from the axial- vector meson Ki. The 
same situation occurs for p and 61 (1235). In Ref. [16], where QCD sum rules for were 
studied, it was argued that one can either explicitly include the contribution of bi in the 
hadronic parametrisation of this "mixed-parity" sum rule and use a suitably large value of 
the continuum threshold Sq ~ 2.1GeV^ or include it in the continuum and use a smaller 
value So ~ 1.0 GeV^. Both procedures yield a stable sum rule and /^(IGeV) ^ 160 MeV. 
For K*, however, the mixed-parity sum rule without Ki does not display a stable plateau in 
M^, which means one has to include the contribution of Ki explicitly.''^ There are actually 
two strange axial-vector mesons, Xi(1270) and fCi(1400), which are usually interpreted as 
mixture of a ^Pi state, the Ka, and a ^Pi state, the Kf, [17, 18]: 

Xi(1270) = KaCOsOK - KbSinOK, 

Xi(1400) = KaSinOK + KbCosOK- 

The results of Refs. [17, 18] indicate that the system is close to ideal mixing, i.e. 6k ~ 45°. 
To the accuracy needed in our sum rules it is then sufficient to replace the two resonances by 
one effective one with the mass = 1.34 GeV [18]. The mixed-parity sum rule obtained 
from the correlation function (7) now reads: 

(/^) v<*/^V (/ij v<i/^^ = 

So So 
m? 



msjss) r /_19 2 

\ 3M^ n \ 9 3 
^Ki also contributes to the sum rule for Jk, but can be safely absorbed into the continuum, as niKi ^ rriK- 
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j2 / ]\/r2 

12M2 ^ TT " ' r M2 V 6 ' V M2 y M2 So V «o 



- 3jgi '".(^'-fG.) - 1^ ((4,)^ + (S.f ) . (10) 

The value of f^^ itself can be obtained from a "pure-parity" sum rule which can be extracted 
from the correlation function (7) leaving one index uncontracted; explicit expressions are 
given in Ref. [14]. As is well-known, the diadvantagc of pure-parity sum rules is that they 
come with a higher mass-dimension which increases the dependence of the result on the 
continuum model, which is however acceptable as long as as we only need an estimate of f^^ 
for use in Eq. (10). The pure-parity sum rule for f^^, including SU(3)-breaking corrections, 
reads [14]: 

so 



so 



Again all scale-dependent quantities are evaluated ji — 1 GeV. The results for /^^ are 
shown in Fig. 2(a), from which we conclude 

(1 GeV) = (0.185 ± 0.010) GeV. (12) 

This result is shghtly smaller than, but still in agreement with, the one obtained in Ref. [19]. 
Using (12) as input in (10), we obtain the values for /^(l GeV) shown in Fig. 2(b), yielding 

/^(l GeV) = (0.185 ± 0.010) GeV. (13) 

The value of at different scales can be obtained from the leading-order renormalisation- 
group improved relation 



4/(3/3o) 

/K(/^) = /i(lGeV) ' 



(^s{iA Y 

as{l GeV) ) 

The result (13) has to be compared with (0.170 ± 0.010) GeV quoted in Ref. [9]. The main 
difference is that the sum rule (10) includes, in addition to the new terms in asfns{ss) and 
G^), in particular the contribution of f^^, which allows one to obtain a stable plateau 
in M^. Once we know the value of f^, we can now determine the continuum threshold to 
be used when Ki is included in the continuum, which will be relevant for the determination 
of a^{K*). Fig. 3 shows that sq = (1.3 ± 0.1) GeV^ is the appropriate value to use if the 
contribution of Ki is not made explicit. 
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Figure 1: (a) fx as function of the Borel parameter for sq = 1.1 GeV^. Solid line: central 
values of input parameters, dashed lines: variation of fx within the allowed range of input 
parameters, (b): the same for /|. with sq — 1.7 GeV^. 
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Figure 2: (a) /^^(IGeV) from Eq. (11) as function of the Borel parameter for sq = 
(2.7, 2.9, 3.1) GeV^ (sohd lines from bottom to top), (b) f^{l GeV) from Eq. (10), using f^^ 
as input. Solid lines: f^^ = 0.180 GcV and, from bottom to top, sq = (2.3,2.45,2.6) GcV^. 
Dashed lines: fj/^_^ = 0.170 GeV (top) and /^^ = 0.190 GcV (bottom). Note that the optimum 
So for f^^ is larger than for in agreement with the resonance structure in the 1"^ and 1~ 
channels. 
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Figure 3: f^{lGeV) from Eq. (10), setting /^^ 
(1.2,1.3,1.4) GeVl 



= 0. From bottom to top: Sq — 
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given in (13) also compares favourably with the average result from lattice calculations 



/^(2GeV) 



= 0.76 ± 0.03 ^ /^(l GeV) = (0.178 ± 0.005) GeV. 



K 



We are now in a position to analyse the sum rules for ai. The sum rule for ai{K) agrees 
with the one obtained in Ref. [6], that for al{K*) is new. We find 

ai(i^)/ie-?^/^' = SRi,+, al{K*){flfe-^KJM^ = SRi,_, 

so 

with SRi ± = A / ds e-/^' (UZ^ 



+ 



18M* ^ TT ' 



5 ms{ss) 



3 M2 



20 a. 



as 

M2 ^ V 



Ei 



So \ , /M^ 



124 8 



27 ^ 9 



1 - 7e + In — + e-'°/^ + Ei 



So 



'^0 

M2 



^ - ms{qq) 5 m,(s(7^G's) 807ra, . ^ _ 

27 TT M2 9 CI i''-^ VW^/ \**/ J ■ 



(14) 



M2 9 81M4 

The sum rules are plotted in Fig. 4; the dominant terms are those in the quark and mixed 
condensates. As the perturbative contribution is ~ O(m^) and hence small, the sum rules 
are not stable in M^, so that we choose to evaluate the sum rules using the optimum values of 
So as determined from the calculation of the decay constants. As expected, the dependence 
of ai on SU(3)-breaking parameters, in particular the precise value of rUs, is much stronger 
than for the decay constants fx- From Fig. 4 we read off the following results for ai. 



ai{K, 1 GeV) = 0.050 ± 0.025, al{K*, 1 GeV) = 0.025 ± 0.015. 



(15) 



The value for ai{K, 1 GeV) agrees with the one obtained in Ref. [6], although our uncertainty 
is shghtly larger. a\{K*) is smaller than ai{K), which follows from the fact that the right- 
hand sides of the sum rules (14) are essentially the same, except for the values of so, and one 
term which gives a positive contribution to ai(K), but a negative to al(K*). Since /|- > fx, 
and the sensitivity of the sum rule on so is small, one clearly expects a\{K*) < ai{K). 
Both results, however, markedly disagree with those obtained in Refs. [4, 5]. In Sec. 1 we 
have already mentioned the reasons for this discrepancy: in Ref. [5] a different set of sum 
rules, nondiagonal sum rules and a chirally odd correlation function, were used which exhibit 
large cancellations between the dominant terms. As discussed in Ref. [6, 7], these sum rules 
are numerically not reliable. Eq. (14) is free of such cancellations and hence expected to be 
more rehable. On the other hand, ai{K) and a\{K*) in (15) are also smaller than the orginal 
results of Chernyak and Zhitnitsky [4]. This is due to the fact that in the corresponding 
nondiagonal sum rule Eq. (6.27) in Ref. [4], which contains no radiative corrections, the 
perturbative term has the wrong sign. Correcting the sign, and using the standard values of 
input parameters from Tab. 1, we obtain the results shown in Fig. 5, which are remarkably 
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Figure 4: (a) ai{K) as function of the Borel parameter for sq — 1.1 GeV^ and = 1 GeV. 
Solid line: central values of input parameters, dashed lines: variation of ai{K) within the 
allowed range of input parameters, (b): Same for al{K*) with sq = 1.7 GeV^. 




Figure 5: ai{K) from the (tree- level) nondiagonal sum rule (6.27) in Ref. [4] after the 
correction of a sign-mistake in the perturbative contribution, sq — (1.0, 1.2, 1.4) GeV^ (from 
top to bottom). 

ai{K*) 




Figure 6: a^{K*, 1 GeV) as function of the Borel parameter for = 1 GeV. Solid hue: 
central values of input parameters, dashed lines: variation of ai{K*) within the allowed 
range of input parameters and Sq — (1.3 ± 0.1) GeV^. 
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close to our result (15) from the diagonal sum rule. Unfortunately, once radiative corrections 
are included in the nondiagonal sum rule, the agreement with (15) is lost and one is back to 
the results obtained in Ref. [5] with negative a^. 

Let us finally turn to ai{K*). As with the decay constant, the sum rule obtained from 
the correlation function (7) contains contributions from the Ki, which in principle need to be 
subtracted. Based on the experience with the sum rule for f^, however, we decide to include 
these contributions in the continuum and use the continuum threshold Sq — (1-3 ±0.1) GeV^ 
determined from Fig. 3. The sum rule for ai{K*) reads 



ml 

QM^V ^ V4 ^ M^J so\so J 



5 msjss) ( rnj_ 

3 M2 1 M2 TT 



(16) 



and is plotted in Fig. 6. We find 

a^(ir*, 1 GeV) = 0.04 ±0.03. (17) 



3 Summary and Conclusions 

In this letter we have calculated the decay constants fx and /|^"'" of the K and K* meson, re- 
spectively, and the first Gegenbauer moments of the leading-twist DAs, ai{K) and a}\^'^{K*). 
We find values for fx and f\ in agreement with experiment and 

/^(l GeV) - (0.185 ± 0.010) GeV, ai{K, 1 GeV) = 0.050 ± 0.025, 

al(K*, 1 GeV) = 0.025 ± 0.015, ai(K*, 1 GeV) = 0.04 ± 0.03. 

The value for agrees, within uncertainties, with that found from lattice calculations [8]. 
The value of ai{K) agrees with that found in Ref. [6], using the same method. ai{K) 
and a\''^{K*) disagree with the negative values found in Ref. [5] which is due to numerical 
instabilities of the nondiagonal sum rules used in that paper. ai{K) and a\{K*) are also 
smaller, by roughly a factor two, than the results obtained in Ref. [7] from exact operator 
relations between ai{K) and a\{K*) and quark-quark-gluon matrix elements. Whereas it 
seems rather unlikely that an increase in the accuracy of the sum rules for oi, Eqs. (14), 
(16), by including corrections in mlas and mgasiqcrgGq) will change the results for oi by 
a factor two, the situation may be different with the sum rules for the quark-quark-gluon 
matrix elements employed in Ref. [7], which only contain terms in ~ rUs- We plan to come 
back to this question in a separate publication. 



10 



Acknowledgements 



R.Z. is supported in part by the EU-RTN Programme, Contract No. HPEN-CT-2002-00311, 
"EURIDICE". 

A Some Details of the Calculation 

In calculating the contribution of the strange quark condensate, one has to include the first 
nonlocal term in the expansion of the quark condensate: 



where i, j are colour, a, (3 spinor indices, and D is the number of dimensions; Q is an auxiliary 
momentum. It is the second term in (A.l) that causes a slight complication in the calculation 
of radiative corrections in the form of finite counter terms. Their origin is twofold: first, the 
factor 1/D induces 0{e) contributions at tree-level which cause finite counter terms upon 
renormalisation. Second, if the derivative in yields a term 7^ in the trace, the contraction 
over K can also yield finite terms in the counter term, which indeed happens for the vertex 
correction diagrams. It appears that these finite counter terms have been missed by the 
authors of Rcf. [6], for we reproduce the results in their appendix using (A.l) and dropping 
just the divergent terms. 
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